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The nature of the dark matter remains a mystery. The possibility of an unstable dark matter 
particle decaying to invisible daughter particles has been explored many times in the past few 
decades. Meanwhile, weak gravitational lensing shear has gained a lot of attention as probe of 
dark energy, though it was previously considered a dark matter probe. Weak lensing is a useful 
tool for constraining the stability of the dark matter. In the coming decade a number of large 
galaxy imaging surveys will be undertaken and will measure the statistics of cosmological weak 
lensing with unprecedented precision. Weak lensing statistics are sensitive to unstable dark matter 
in at least two ways. Dark matter decays alter the matter power spectrum and change the angular 
diameter distance-redshift relation. We show how measurements of weak lensing shear correlations 
may provide the most restrictive, model-independent constraints on the lifetime of unstable dark 
matter. Our results rely on assumptions regarding nonlinear evolution of density fluctuations in 
scenarios of unstable dark matter and one of our aims is to stimulate interest in theoretical work on 
nonlinear structure growth in unstable dark matter models. 
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I. INTRODUCTION 

A preponderance of evidence supports a picture in 
which ^ 5/6 of the mass density of the Universe resides 
nonbaryonic dark matter (reviews include Refs. [IHll)- 
The prevailing hypothesis is that the dark matter is an 
as yet undetected particle that survives as a relic from the 
hot, early Universe. An effort to identify the dark matter 
now proceeds on many fronts and the dark matter is cur- 
rently constrained by direct searches (e.g., Refs. 
indirect searches (e.g., Refs. |^, [To|). and astronomical 
observations (e.g., Refs. |ll| - [l4| ). In this paper, we ex- 
plore the possibility of constraining invisible decays of the 
dark matter particle using forthcoming statistical mea- 
surements of weak gravitational lensing. 

Limits on unstable dark matter have been considered 
by numerous authors in the recent literature [l5l - [25j . Ra- 
diative decays are very strictly limited, with the best 
constraints yielding lifetime bounds of tddmJc^ 10''iJQ"^ 
[26l - l30| . Assuming decays to Standard Model neutri- 
nos, the least detectable Standard Model particles, places 
mass-dependent limits as restrictive as tddm 
for dark matter particle masses near 10^ GeV [2^. [sij. 
though constraints are strongly mass dependent. Cos- 
mological tests provide an opportunity to constrain the 
stability of the dark matter independent of particle mass 
and the interactions of the decay products, and current 
cosmological limits on invisible dark matter decays imply 
that the lifetime for decays to relatively light products is 

We consider the possibility of improving model- 
independent constraints on the dark matter particle life- 
time using forthcoming weak lensing data. Such inde- 
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pendent constraints would be most relevant to models of 
light dark matter (masses ^10 MeV to evade neutrino 
constraints) or dark matter sufficiently seq uestered from 
the Standard Model (e.g., Refs. [H [ll, [sJil]) and may 
help improve or complement constraints of asymmetric 
dark matter models 31, 36, Is^l- For concreteness, we 
consider constraints on dark matter lifetime in a bench- 
mark model of a cold dark matter particle that undergoes 
two-body decay to light daughter particles with a lifetime 
tuned to be exceptionally large. 

The primary constraint comes from scale- and redshift- 
dependence of the cosmological gravitational lensing 
power spectrum at z ~ — 3 after normalizing the power 
spectrum of density fluctuations at z w 1100 via cos- 
mic microwave background (CMB) measurements. A 
full exploration of possible constraints is difficult be- 
cause the nonlinear evolution of structure in the Uni- 
verse in such models has not been extensively studied. 
In our most conservative forecasts for what may be pos- 
sible with forthcoming instruments, we find that utiliz- 
ing only scales on which linear perturbative evolution 
should be useful {£ < 300) and taking weak prior con- 
straints on other cosmological parameters, forthcoming 
surveys may produce dark matter lifetime constraints 
that are, at minimum, competitive with contemporary, 
model-independent constraints. More aggressive priors 
expected from Planck CMB measurements improve upon 
these constraints by roughly a factor of two. We ar- 
gue that utilizing the information that may be available 
from nonlinear evolution should improve upon these con- 
straints by an order of magnitude. Our most aggressive 
forecasts, using information extending to scales £ < 3000 
and taking Planck CMB prior constraints on other cos- 
mological parameters, suggest that constraints as strong 
as tddm^ a few x IO^Hq^ may be possible with a sur- 
vey covering a large fraction of the sky, such as that 
planned by the Large Synoptic Survey Telescope (LSST) 
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[38| or Euclid [33|- Achieving rehable constraints from 
such measurements requires an understanding of nonlin- 
ear evolution in unstable dark matter models. The pos- 
sibility of stringent constraints on decaying dark matter 
from astronomical imaging surveys should be strong mo- 
tivation to study the nonlinear evolution of cosmic struc- 
ture formation in such scenarios (see also Refs. 

We continue our manuscript with a discussion of weak 
lensing observables in § [iTl We describe our benchmark 
model for unstable dark matter and the evolution of cos- 
mological perturbations in such a model in § IIIII We 
discuss nonlinear structure evolution and the two meth- 
ods we use to estimate nonlinear evolution in § IIVI We 
describe our methods for forecasting constraints on un- 
stable dark matter in § |Vl This section also includes 
a summary of our fiducial cosmological model and our 
assumptions regarding prior constraints. We illustrate 
the effects of unstable dark matter on lensing observables 
and present our forecast limits on dark matter lifetimes 
in § IVIl Finally, we summarize and discuss avenues for 
future work in § IVIII 



II. 



WEAK GRAVITATIONAL LENSING 
OBSERVABLES 



We explore the utility of weak gravitational lensing 
measurements for constraining the stability of the dark 
matter. Our most robust forecasts derive from considera- 
tions of possible weak lensing measurements restricted to 
scales where linear perturbative evolution of the metric 
potentials remains useful. In this manner, our paper is 
very similar in spirit to that of Schmidt (40| , who studied 
constraints on modified gravity from weak lensing statis- 
tics restricted to linear scales. However, we attempt to es- 
timate possible improvements to the constraining power 
of weak lensing observables, provided that nonlinear evo- 
lution can be modeled robustly. 

We consider the set of observables that may be avail- 
able from ongoing and forthcoming large-scale galaxy 
imaging surveys to be the auto- and cross-spectra of lens- 
ing convergence from sets of galaxies in ^Vtom redshift 
bins. The iVTOM(^TOM + l)/2 distinct convergence spec- 
tra are 



where i and j label the redshift bins of the source galax- 
ies. We take iVxoM = 5 and consider evenly-spaced bins 
in redshift from a minimum redshift of z = to a max- 
imum redshift of z = 3. In agreement with the study of 
Ma et al. [41jj , we find that finer binning is not required to 
maximize the constraining power of such surveys. Weak 
lensing as a cosmological probe has been discussed at 
length in numerous papers (a recent review is Ref. p^). 
We give a brief description of our methods below, which 
are based on the conventions and notation in Ref. (isj 
(to which we refer the interested reader for details). 



In Eq. [1] H{z) is the Hubble expansion rate. Da is the 
comoving angular diameter distance, and Pit-<i)(fc,z) is 
the power spectrum of Newtonian gauge scalar potentials 
^ — $ at wavenumber k and redshift z. In the follow- 
ing section, we describe our use of the publicly-available 
CMBFAST code to calculate P*_$(fc,z), in which case it 
will be more natural to work in the synchronous gauge. 
Transforming between coordinate systems can be accom- 
plished straightforwardly by following, for example, the 
methods described in Ref. [HI which we do not repeat 
here. 

The Wi are the so-called lensing weight functions for 
source galaxies in redshift bin i. In practice, the galaxies 
will be binned by photometric redshift, so that the bins 
will have nontrivial overlap in true redshift (see Ref. [ill 
for a detailed discussion). Defining the true redshift dis- 
tribution of source galaxies in the ith photometric red- 
shift bin as drii/dz, the window functions are 



W,{z)^Da / dz 



j Da{z,z') duj 
Da{z') dz' 



(2) 



where Da{z^ z') is the angular diameter distance between 
redshift z and z'. 

We model the increased uncertainty induced by uti- 
lizing photometric galaxy redshifts with the probability 
function of assigning an individual source galaxy photo- 
metric redshift Zp given a true redshift z, P{zp\z). In this 
notation, the true redshift distribution of sources in the 
ith photometric redshift bin is 



dni{z) 
dz 



dZr 



dn{z) 
dz 



P(zp|z) 



Here we take the true redshift distribution to be 



dn{z) 
dz 



4z^ 



V2^ 



exp[-(z/zo)^ 



(3) 



(4) 



with zq ~ 0.92, so that the median survey redshift to 
Zmed = 1, and n as the total density of source galaxies 
per unit solid angle p5l - l47l |. We assume that uncertain 
photometric redshifts can be approximated by taking 



P(zp\z) = 



1 



: exp 



(zp - z)2 



2a2 



(5) 



where az{z) — 0.05(1 -I- z) [4l|. Complexity in photo- 
metric redshift distributions is something that must be 
overcome to bring weak lensing constraints on cosmology 
to fruition (e.g., Ref. [ii,!!^). 

Observed convergence power spectra P^-'{i), contain 
both signal and shot noise. 



p:^{i) = p^ + nA,(rf') 



(6) 



where (7^) is the noise from intrinsic ellipticities of source 
galaxies, and Ui is the surface density of galaxies in the 
ith tomographic bin. We follow the recent convention and 



3 



set ^ (7^) = 0.2, subsuming additional errors on galaxy 
shape measurements into an effective mean number den- 
sity of galaxies, n. Assessments of intrinsic shape noise 
per galaxy may be found in, for example [s^ Isol . [Slj . 
Assuming Gaussianity of the lensing field, the covariance 
between observables Pl^ and P^^ is 



pik p]l pil pjH 



(7) 



where the i and j map to the observable index A, and 
k and 1 map to B such that Cab is a square covariance 
matrix with -/Vtom(A'tom + 1)/2 rows and columns. We 
assume Gaussianity throughout this work and consider 
only multipoles £ < 3000 at which point the Gaussian as- 
sumption and several weak lensing approximations break 

down lii-ii. 



III. PERTURBATION THEORY WITH 
UNSTABLE DARK MATTER 

Our aim is to predict the power spectrum of weak grav- 
itational lensing convergence in unstable dark matter sce- 
narios. To do so, we must compute the modifications 
to the metric potentials in unstable dark matter scenar- 
ios [see Eq. ([T])]. We explore a restricted set of models 
in which a massive parent dark matter particle decays 
into a significantly lighter pair of daughter particles. For 
the sake of specificity, we adopt a decaying dark mat- 
ter (DDM) scenario in which massive majorana parent 
particles decay into relativistic daughter (RD) particles 
via two-body decay and use this scenario to benchmark 
observational constraints. In such a scenario, the life- 
time of the unstable dark matter particle lifetime (r) is 
the only nonstandard free parameter. One could assume 
decay to a combination of heavy and light daughter par- 
ticles in which case the mass differences are important 
additional parameters that establish the recoil velocities 
of the decay product particles (as explored in detail in 
Refs. recently). 

The distribution functions of DDM (/ddm) and 
RD (/rd) evolve according to the coupled Einstein- 
Boltzmann equations. In particular (e 



coupiea Jim 

■g-, Ullii), 



DDM 



dr 

d/RD 



tddmeddm ' 



/DDM 



tddm 



/DDM 



dr 



tddmeddm ■ 



/] 



1 



DDM 



'''DDM 



-/l 



DDM, 



(8) 
(9) 



where r is the conformal time and tddm, ^ddm, and 
TTT-DDM are the lifetime, energy, and mass of decaying 
dark matter. Following established procedure, we express 
the distribution function of species X as a zeroth-order 
distribution plus a perturbation. 



fxix,q,T) /^(g,T)[l + *x(:r,<f,r)] 



(10) 



The evolution of the mean energy density for DDM and 
its RD particles follow from the zeroth-order integrals of 



Eq. dH) and Eq. 



PDDM + 3-ffpDDM 
PRD + 4i7/7RD 



— FpDDM 
FpDDM 



(11) 

(12) 



Here and throughout, we designate y as the time deriva- 
tive of y, and we denote the decay rate as F = 1/tddm. 
In the limit of a massive DDM particle, evolution of the 
comoving density pbdmo-'^ approaches cx exp(— ^/tddm)- 
The collision term describing the DDM decays is pro- 
portional to /ddm> rendering the equations describing 
the evolution of DDM perturbations identical to those 
of standard, stable cold dark matter at the lowest order 
in perturbation theory. The perturbation equations de- 
scribing the daughter particles are less trivial. Following 
Refs. 15, 44], we expand the perturbation equations for 
RD particles in a series of Legendre polynomials Pi{x), 
yielding 



-FRD(A:,n,T) 



Jdqq^fl^^iq, r)*RD 



Jdqq^f^j,{q,T) 

00 

= ^(-i)'(2/ + l)FniyAk, T)Pi{k- n), (13) 



/=o 



where F^j^iik^r) are the harmonic expansion coeffi- 
cients. The orthonormality of Legendre polynomials al- 
lows the evolution equations to be written as 



^RD = 



„^ X , Pddm,. 

^prdJ H (.(JDDM 



PRD 

3 , 2 / '^RD N PDDM , 

7RD = (— ; ctrdJ 1 

4 PRD 



7RD 



<^rd) (14a) 
(14b) 



CRD 
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(261; 



RD 



h + Qf] 



-kF, 



RD,3j 

Pddm 

PRD 



-O-RD 



(14c) 



F 



RD 



21 + 1' 



m 



RD,/-1 



(/ + 1)-Frd,/+i 
Pddm 



Prd 



-Frd./, I > 3 (14d) 



at first order, where (5rd = Frd,i, 6'rd = 3/4A;Frd,i, 
ctrd = ^RD,2, and h is the scalar trace of the metric 
perturbation, all in well-established notation. 

We have modified the publicly- available CMBFAST code 
of Seljak and Zaldarriaga [57] to compute the poten- 
tial power spectra. As we noted in § [TTl we quote the 
perturbation equations explicitly in synchronous gauge 
simply because CMBFAST is written in terms of the syn- 
chronous gauge perturbations. Gauge transformations 
can be made straightforwardly [43 |. 

The growth of perturbation is affected by the change 
of energy density among the relativistic and nonrelativis- 
tic components. From Eq. (ITT|) and Eq. we can see 
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that in the decaying dark matter scenario the comov- 
ing dark matter density decreases exponentiahy, and ah 
of this decrement is transferred into relativistic energy 
density. Consequently, perturbation growth exhibits a 
scale-dependent suppression relative to stable dark mat- 
ter, where the relevant scale is the horizon size at the 
epoch of decay. This late-time suppression of structure 
growth in large part provides the necessary leverage for 
weak lensing constraints on unstable dark matter. For 
daughter particles, the additional Pddm I Prd terms have 
an impact on scales greater than the horizon at the time 
of decay [11] . 

As we will discuss below, some of the constraining 
power of weak gravitational lensing comes from obser- 
vations made on scales where linear perturbation theory 
is no longer adequate (e.g., [M El [H, [H-li^l). The 
constraints we forecast in the following sections that are 
based on linear scales only are robust and interesting in 
and of themselves. However, the utility of weak lensing 
is greatly increased if scales modified by nonlinearity can 
also be exploited for cosmological constraints 43] , so we 
explore multiple proposed nonlinear corrections to linear 
evolution in the following section. 



IV. NONLINEAR EVOLUTION 

In the standard application of weak lensing to con- 
strain dark energy, most of the constraining power comes 
from scales on which linear evolution of cosmological per- 
turbations is no longer valid. Accounting for nonlinear 
evolution enables a larger range of multipoles to be used, 
and nonlinear evolution greatly enhances signal-to-noise 
of weak lensing measurements at multipoles £^ 300. To 
estimate the constraints that may be anticipated from 
a full, nonlinear treatment of DDM, we explore nonlin- 
ear corrections to linear evolution using both the method 
of Smith et al. [6l|, and a halo model-inspired method 
by Peter [l^. Smith et al. [6l| provide an empirical fit 
for nonlinear power given a linear power spectrum. We 
utilize the fit of Smith et al. [6l| directly as one of our 
nonlinear structure models. This is not entirely unrea- 
sonable, because lifetimes of interest are far larger than a 
Hubble time, so little decay occurs relative to a standard 
cosmological model. We implement the method of Peter 
[ll] using the halo model as follows. 

The halo model (see Ref. jg^l for a review) is based 
on the assumption that all matter resides within dark 
matter halos. The matter power spectrum is given by 
the sum of two terms, 



and 



(15) 



where 



P2H{k) 



1 

pIi 



dm m— — X{k\m) bh{m) 
dm 



(17) 

In the foregoing equations, pm is the mean matter den- 
sity of the universe, m is halo mass, X{k\m) is the Fourier 
transform of the Navarro, Frenk, and White (NFW, 
Ref. 63]) density profile for a halo of mass to, P''"(/c) 
is the linear matter power spectrum, and bh{m) is the 
halo bias function. The one-halo term Pijj(/e), describes 
correlations among mass elements within a common halo 
while the two- halo term P2H{k), is due to correlations 
among mass elements in distinct halos. 

To estimate the impact of decaying dark matter on 
matter clustering, we follow the approach denoted as 
Case 1 by Ref. [la] to describe modifications to the halo 
mass function, halo bias, and internal halo structure. We 
then incorporate these modifications into the halo model 
formulae of Eq. ([T6l) and Eq. ([T7|) to compute lensing 
power spectra. This model is based upon the assump- 
tion that halos at early times are very much like their 
counterparts in models of stable, cold dark matter (be- 
cause little decay will occur in any viable model) and 
that modifications to halo structure can be described by 
the conservation of adiabatic invariants describing dark 
matter particle orbits. 

Consider a population of dark matter halos that 
formed prior to any significant dark matter decays such 
that halos at any time t <C H^^ can be modeled as stan- 
dard, CDM halos. These halos then lose mass as their 
constituent dark matter particles decay. If the decay life- 
time is much larger than the halo dynamical timescale (as 
it will always be in cases of interest because dynamical 
times are Tdyn ^ Q.IHq^ and viable regions of parameter 
space are tddm ~> Hq^)- then the halo gravitational po- 
tential changes adiabatically. Exploiting the adiabatic in- 
variance of angular momentum for particles on nearly cir- 
cular orbits, establishes a prediction for the relationship 
between the initial and final matter distribution within 
a dark matter halo, 



(18) 



Pi^(fc) = ^ / dmm'^\^{k\m), (16) 
Pli J dm 



where Mi{r) is the mass enclosed within radius r in the 
initial, early-time halo, Mf(r) is the corresponding quan- 
tity describing the contemporary, late-time halo, and 
and rf are the initial and final radii of a particle shell, as- 
suming that mass shells never cross and particles move in 
circular orbits. Eq. (fT8|) is the basic relation of the stan- 
dard, adiabatic contraction model for predicting mod- 
ifications of halo structure due to coUisional processes 

[iiii]. 

For unstable dark matter, with a lifetime Tddm, a frac- 
tion f{Tddm, z) of unstable dark matter particles will have 
decayed by redshift z. According to the adiabatic con- 
traction model, the mass enclosed in rj will be 



Mfirf) - (1 - fiTddm,z))M,{n) 



(19) 
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Inserting Eq. ^ into Eq. (|T8 
the initial and final radii is 



the relationship between 



(20) 



If we assume that the initial dark matter halos can be well 
described by NEW profiles, the final mass distribution 
will be 



1 dM 



/ 



47rr^ drf 
(1 - /)2 dM, 



dvi 
(1 " JYPs 



(1 - f)rf 



1 



(1 - fYf 



(21) 
(22) 
(23) 



We model the initial mass function drii/dm and halo 
bias hh (m) using the relations of Ref. [HI] ■ This choice is 
made for convenience because in models with stable dark 
matter, it satisfies the necessary conditions that the halo 
model integrals contain all mass and that the clustering 
of dark matter is unbiased with respect to itself. Some 
definitions of halo virial radii will be altered by decays. 
In order to ensure that all mass remains accounted for, 
we define halos as the mass within virial radii fixed to a 
definition of 200 times the average density of the Universe 
in the absence of decays. Thus, virial radii are fixed to 
be the same as they would be in standard CDM, but halo 
masses are smaller by a factor of 1 — /(tddm,^:)- This 
definition preserves the convenient properties of the bias 
and mass relations in Ref. [g^] and is identical to their 
halo definition in the absence of dark matter decays. 

The new halo mass function at mass Mf is 



dn/(M/[M,],z) _ dn,(M,,z) 



dM/ 



dMi 



dM; 



dM/ 



and 



(24) 



(25) 



where the initial and final masses are related via Eq. (1191) . 
In other words, we assume the abundance and clustering 
to follow the abundance and clustering laws for halos 
of stable dark matter of the corresponding masses. No- 
tice that the abundance of halos of a given contemporary 
mass Mf is reduced compared to that in a stable dark 
matter model because the final mass refiects the mass 
loss due to decays and more massive halos are intrinsi- 
cally rare. Likewise, halos of final mass Mf are more 
strongly clustered than their counterparts in stable dark 
matter scenarios because halo bias is an increasing func- 
tion of mass (see Ref. |68j for the basic theory of the 
mass function and bias). The halo density profiles also 
become shallower as Ts increases and ps decreases when 
the decay-induced modifications to halo profiles are ac- 
counted for. 
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FIG. 1: Fractional difference of shear lensing auto-spectra 
with and without halo mass loss due to dark matter decay. In 
this panel, we choose a dark matter lifetime of r = 500 Gyr 
for illustrative purposes. From top to bottom, lines are for 
the auto-spectra of sources in the first tomographic redshift 
bin (0 < 2:p < 0.6) to the fifth tomographic redshift bin 
(2.4 < Zp < 3.0). In this panel, both spectra are computed 
using a halo model. For comparison, note that the differences 
between the halo model and the Smith et al. [gJI relation 
can be as much as 20% with the maximum discrepancy near 

e ~ 1000. 



The reduction in the number of halos and the shal- 
lowing of halo profiles reduces lensing power compared 
to a halo model with no accounting for mass loss. Eig- 
ure [T] shows a comparison between halo model calcula- 
tions of lensing power spectra including and neglecting 
halo mass loss. The greatest changes are at relatively 
high £ {£ ^ 300) and are due to the concentration change 
which alters the one-halo term [Eq. p^ ] . The shift in the 
mass function and halo bias cause the slight reduction in 
the two-halo term and power at lower £. As we show in 
§ IVIi this additional reduction in power is a distinctive 
feature that leads to slightly more restrictive bounds on 
DDM lifetimes. 

We emphasize that neither of these approaches have 
been calibrated in detail on simulations of structure 
formation in DDM cosmological models. However, we 
demonstrate that the region of parameter space relevant 
to forthcoming constraints has tddm ^ (see also 

[T^-i^l)- This means that little of the DDM wiU have de- 
cayed prior to the present epoch and the boost in signal- 
to-noise should be something close to that afforded by the 
nonlinear treatment of standard, stable dark matter. In 
actuality, only a detailed numerical treatment can answer 
these questions definitively (as Refs. [3 13 have recently 
argued). It is our hope that this proof-of-concept paper 
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will motivate pursuit of large-scale simulations of DDM 
similar to those being carried out in support of the dark 
energy constraint program (e.g., |69l471j ). 



V. FORECASTING METHODS 

The Fisher Information Matrix provides a simple es- 
timate of the parameter covariance given data of spec- 
ified quality. The Fisher matrix has been utilized in 
numerous, similar contexts in the cosmology literature 
[43I lio, [zl-lzl] , so we give only a brief review of im- 
portant results and the caveats in our particular applica- 
tion. 

The Fisher matrix of observables in Eq. ([T]), subject to 
covariance as in Eq. ([7]), can be written as 



E 



A,B 



dpi 



AB- 



Ft 



(26) 

where the indices A and B run over all iVTOM(-^TOM + 
1)1 2 spectra and cross spectra, the pi are the parameters 
of the model, fsky is the fraction of the sky imaged by the 

— 1/2 

experiment, and imin = ^f^/^y is the smallest multipole 
constrained by the experiment. F,^ is a prior Fisher ma- 
trix incorporating previous knowledge of viable regions 
of parameter space. We set imax = 3000 in our most 
ambitious forecasts. On smaller scales (higher £), vari- 
ous assumptions such as the Gaussianity of the lensing 
field, break down [43, 52- 56, 79]. To be conservative, we 
explore modest priors to each parameter independently, 
so that F^ = Sij/{af)'^, where erf is the la prior on 
parameter The forecast, la, mar ginalized constraint 
on parameter pi is a{pi) = \/[F~^ii. 

Other than the DDM lifetime tddm, we vary six cos- 
mological parameters that we expect to modify weak 
lensing power spectra at significant levels and to exhibit 
partial degeneracy with tddm- We construct our fore- 
casts for DDM lifetime bounds after marginalizing over 
the remaining parameters. Our six additional param- 
eters and their fiducial values (in parentheses) are the 
dark energy density (0.74), the present-day dark mat- 
ter density, wdm = f^DM^^ (0.11), the baryon density 
ujb ~ ribh^ (0.023), tilt parameter ns(0.963), the natural 
logarithm of the primordial curvature perturbation nor- 
malization In(A^) (—19.94), and the sum of the neutrino 
masses J^i "^i^i (0.05 eV). This implies a small-scale, low- 
redshift power spectrum normalization of erg ~ 0.82. The 
optical depth to reionization has a negligible effect on the 
lensing spectra on scales of interest, so we do not vary 
it in our analysis. We adopt as our null hypothesis a 
stable dark matter particle with F = 0. Note that when 
r 7^ 0, there is a higher density of dark matter in the 
past than would be inferred for a stable dark matter par- 
ticle because we choose our parameter set to describe the 
contemporary matter density. 

We take priors on our cosmological parameters of 



cr(a;„0 = 0.007, cr(wh) = 1.2 x lO'^, a{\nAj^) = 0.1, 
a{ns) — 0.015, and cr(f2A) — 0.03. We assume no pri- 
ors on DDM lifetime or neutrino mass. Our fiducial 
model is motivated by the WMAP seven-year result and 
our priors represent marginalized uncertainties on these 
parameters based on the WMAP seven-year data [so! ]. 
These priors are conservative and allow for weaker con- 
straints on DDM than would be expected from future 
data, where stronger priors may be available. To esti- 
mate the potential power of lensing constraints on DDM 
when stronger cosmological constraints are available, we 
also explore prior constraints on these parameters at the 
level expected from the Planck mission^ using the entire 
Planck prior Fisher matrix of Ref. [8l|. Of course, using 
published priors from other analyses is not self-consistent 
because these priors were derived in analyses that assume 
stable dark matter, but relevant lifetimes the dark matter 
decays should cause only subtle alterations to the cosmic 
microwave background anisotropy spectrum so this anal- 
ysis should approximate a self-consistent simultaneous 
analysis of all data. 

In some cases, we will estimate nonlinear power spec- 
tra in models with significant neutrino masses. In such 
cases, we follow the empirica,l prescription established in 
previous studies (e.g., Refs. (zJlH,!!^) and take 



where 



/b-l-DM — 



DM 



(27) 

(28a) 
(28b) 



P^'"(fc) is the linear power spectrum of neutrinos, and 
^b+bivil^) the nonlinear power spectrum evaluated for 
baryons and dark matter only. Again, our adoption of 
this prescription may induce errors in our calculation and 
only a large-scale numerical simulation program can test 
this assumption. 

We explore possible constraints from a variety of forth- 
coming data sets. We consider the Dark Energy Survey 
(DES)^ as a near-term imaging survey that could pro- 
vide requisite data for this test. We model DES by tak- 
ing a fractional sky coverage of fskij — 0.12 and with 
n = 15/arcmin^. Second, we consider a comparably nar- 
row, deep imaging survey as may be carried out from a 
space-based platform, such as a Supernova Acceleration 
Probe-like implementation of a Joint Dark Energy Mis- 
sion (JDEM)'^''* or the National Academy of Science's 



^ http://www.esa.int/planck 

^ http : //www . darkenergy survey . org| 

^ http: //universe .nasa.gov/p rogram/probes/jdem I 
^ http://snap.lbl.gov/ 
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Fisher matrix formalism cannot formally be applied to 
parameters near physical limits in their parameter values. 
In such Fisher approach allows for parameter de- 

generacies that extend into the forbidden region of the 
parameter space and should not be permitted on physi- 
cal grounds. This additional degeneracy tends to cause 
under estimates of constraints that may be realized from 
a more detailed analysis. An example of this is neutrino 
mass. The Fisher matrix has been utilized to constrain 
neutrino mass and empirically Fisher matrix projections 
for neutrino mass constraints have been shown to match 
well direct searches of parameter space (e.g., [ttI. [s^ - lsBj ) 
To verify our Fisher matrix results, we have performed 
several direct searches of subspaces of our full parameter 
space (limited by computational cost) to constrain DDM 
lifetime. Generally, we find the marginalized Fisher con- 
straints to be only slightly less constraining than the di- 
rect search results and we present an example of this in 
the following section. The computational cost of a full 
parameter search seems unwarranted given the theoreti- 
cal limitations discussed in § IIIII 



FIG. 2: Relative diflerence in the linear potential power spec- 
tra between DDM and stable dark matter models a,t z — 4 
(top) and z = (bottom). The solid lines show a model with 
tddm = 10^ Gyr while the dashed lines have tddm ~ 300 Gyr. 
The dash-dotted lines show the influence of a non-negligible 
neutrino mass with = 0.3 eV. 



Astronomy and Astrophysics Decadal Survey suggestion 
of a Wide Field Infrared Space Telescope (WFIRST)5. 
We refer to such a survey as a "Deep" survey and model 
it with fsky = 0.05 and fi = 100/arcmin^. Lastly, we 
consider a class of future "Wide" surveys as may be car- 
ried out by the Larg e Synoptic Survey Telescope (LSST)^ 
[ssj or Euclid'' We model these Wide surveys with 
fsky = 0.5 and fi = 50/arcmin^. In all cases, we take 
(7^) = 0.2 and assume particular shape measurement 
errors from each experiment are encapsulated in their 
effective number densities, in accord with recent conven- 
tional practice in this regard. Our results are relatively 
insensitive to number density because shot noise does not 
dominate cosmic variance on the scales we consider for 
any of the experimental parameters we consider. 

The Fisher matrix is valuable because it greatly re- 
duces the computational effort necessary to forecast con- 
straints from forthcoming experimental data. However, 
the Fisher matrix formalism has important drawbacks. 
First, the Fisher matrix only characterizes parameter de- 
generacies locally about the fiducial model. Second, the 



^ 'http : // sites . nationalacademles . org/bpa/BPA_049810| 
^ http : // www. Isst ■ org| 
[http : //sci .esa. int/euclld| 



VI. RESULTS 
A. Weak Lensing Power Spectra 

Weak lensing power spectra are altered by DDM in 
two respects. First, the power spectra for potential and 
density fluctuations are altered in a scale-dependent way. 
At early epochs, when the matter density is higher in 
the DDM models than in standard ACDM, potential 
and density fluctuations are larger because the epoch of 
matter-radiation equality occurs earlier. We have verified 
that our constraints are insensitive to the epoch at which 
we normalize the matter density. At late times, DDM 
decays suppress density and potential fluctuations. We 
show this dependence of potential fluctuations on DDM 
lifetime in Figure [H Notice that models of unstable dark 
matter have greater P*_$(fc) on scales 10^^ h Mpc^^ 
at high redshift, but this power is suppressed on subhori- 
zon (k > 10~^ h Mpc~^) at lower redshifts. The strong 
scale dependence in potential power spectra at scales of 
order k ~ 0.05 h Mpc~^ should be present in convergence 
spectra projected on multipoles £ ~ fc Da{z = I) '-^ 150 
(z = 1 is the median redshift of lensed sources in our 
model surveys). The different redshift dependence of 
DDM, which results in greater suppression of power with 
decreasing redshift, compared to neutrino mass-induced 
power suppression allows the two to be disentangled. 

The observed strength of gravitational lensing also has 
a dependence upon geometry, so differences in angular 
diameter distance may lead to modified lensing power 
spectra. These geometrical differences provide the bulk 
of the information with which lensing can constrain dark 
energy [sl, H^l and the angular diameter distance to the 
last-scattering surface has been used to constrain decay- 
ing dark matter in previous studies [itI [20j . In princi- 
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FIG. 3: Relative difFerences in comoving angular diameter 
distance between models with DDM and stable dark matter 
as a function of scale factor. The solid line represents a model 
with tddm = 10^ Gyr and the dashed line shows a model with 
TDDM = 300 Gyr. 



pie, the distance-redshift relation in decaying dark mat- 
ter models can be mimicked by dark energy with a vari- 
able equation of state [sll . We show in Fig. [3] a compar- 
ison of the angular diameter distance in DDM models. 
Fig. [3] demonstrates the angular diameter distances are 
modified at levels that are small compared to the rela- 
tive potential fiuctuations shown in Fig. As a conse- 
quence, we find that DDM constraints are based mainly 
on the scale-dependent potential fluctuation modifica- 
tions, rather than on the modified distance scale which 
is the primary driver of dark energy constraints. 

These changes manifest in a scale- and redshift- 
dependent shift in observable convergence power spectra. 
We show examples of the shift in the convergence power 
spectra in models of unstable dark matter in Figure SI 
Notice the strong scale dependence on multipoles of a 
few hundred. On smaller scales, the net effect is an over- 
all change to the level of the convergence power, with 
only a weak scale dependence. This is one reason that 
we anticipate that we may be able to utilize the methods 
of Smith et al. [6l| or Peter [l^ to approximate the con- 
straining power of weak lensing surveys when nonlinear 
evolution is included. 



B. Forecast Constraints on DDM Lifetimes 

As a first attempt to estimate the power of weak lens- 
ing to constrain DDM, we examine models in which we 
consider only the linear evolution of potential and density 
fiuctuations. We consider taking the maximum multipole 
that we may observe as high as imax = 3000; however. 
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FIG. 4: Relative modification to linear lensing convergence 
power spectra caused by DDM. The lines show the conver- 
gence power spectra for galaxies in our third redshift bin 
(galaxies with photometric redshifts 1.2 < Zp < 1.8) in mod- 
els with finite tddm compared to a model with stable dark 
matter. The solid line shows the power spectrum residual in 
a model with tddm = 1000 Gyr, and the dashed line shows 
the same power spectrum residual with tddm = 300 Gyr. 



for multipoles larger than £}Z 300 nonlinear effects will 
be very important (also, see Fig. [Ij. Therefore, we 
also consider taking £rnax = 300 so that we consider only 
those scales for which linear perturbative evolution may 
be valid. In this case, our theoretical methods are appli- 
cable and forthcoming weak lensing constraints should 
do at least this well. 

In our most ambitious forecasts, we assume that we 
can use the nonlinear formula of Smith et al. [61.] or the 
halo model to estimate the boost in signal that nonlinear 
evolution may provide for weak lensing constraints on 
DDM. These approaches toward nonlinear corrections are 
not entirely self-consistent, but may serve as an indicator 
of what could be achieved if a numerical simulation effort 
addressed nonlinear evolution in DDM robustly. 

We summarize our primary results for the upper lim- 
its that may be set on the DDM decay rate F, by weak 
lensing measurements in Table U The limits in this table 
have been marginalized over all other cosmological pa- 
rameters, including neutrino masses. We computed the 
results in the upper portion of Table |I] using contempo- 
rary priors on other cosmological parameters. Results 
below the middle dividing line of Table |T] were computed 
with prior constraints on cosmology at levels expected 
from the Planck CMB mission and are labeled with a 
"PP." Different lines in Table |T] show results using differ- 
ent model power spectra. The options are the linearly- 
evolved power spectrum only, results correcting for non- 
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TABLE I: Forecast 68% marginalized limits on dark mat- 
ter decay rates from weak lensing surveys under several as- 
sumptions. The limits are in units T/Hq, where Hq — 
72 km/s/Mpc. Constraints are shown for "Linear" power 
spectra, "Smith et al." nonlinear corrections, "Halo Model" 
nonlinear corrections, and "Modified Halo Model" nonlinear 
corrections that account for mass loss as in Ref. [f^. The 
abbreviation "PP" stands for Planck priors. 



Experiment 



DES Deep Wide 



Linear, imax = 3000 0.07 0.06 0.046 

Linear, l^ax = 300 0.08 0.09 0.057 

Smith et al, £„,ax = 3000 0.03 0.02 0.008 

Smith et al, ln,ax = 300 0.06 0.05 0.029 

Halo Model, i^ax = 3000, 0.03 0.02 0.010 

Modified Halo Model, £,nax = 3000, 0.02 0.02 0.008 

Linear, Imax = 3000, PP 0.03 0.03 0.016 

Linear, Imax = 300, PP 0.06 0.07 0.026 

Smith et al, £,nax = 3000, PP 0.02 0.01 0.006 

Smith et al, imax = 300, PP 0.05 0.05 0.018 

Halo Model, £,^ax = 3000, PP 0.02 0.02 0.007 

Modified Halo Model, £,nax = 3000, PP 0.02 0.01 0.006 



linear evolution using the Smith et al. |6l| formula, non- 
linear power results using the halo model, and nonlinear 
power using the halo model modified to account for the 
loss of mass within halos (following Ref. [l^)- In each 
case, we consider both restricting to linear scales tak- 
ing ^max = 300 and using nonlinear information with 
^max — 3000 to constrain decaying dark matter. 

Constraints that exploit only linear scales are already 
promising. A DES, a Deep JDEM/WFIRST-like survey, 
or a Wide LSST- or Euclid-like survey should constrain 



the DDM lifetime at the level of tddm = T 



-i> 



13m 



12Hq , and 18Hq with contemporary priors on other 
cosmological parameters. These results are already com- 
parable to contemporary, model-independent constraints 
on unstable dark matter ^- 21] and do not require de- 
tailed calibration of nonlinear structure growth or ambi- 
tious priors on other cosmological parameters (In A|j in 
particular) . It seems reasonable then, that weak gravita- 
tional lensing will provide, at minimum, a complemen- 
tary, model-independent technique to constrain DDM 
that is competitive with other, existing techniques. 

If we interpret the other entries in Table U as pos- 
sible limits that may be achieved if the necessary 
nonlinear evolution in models of DDM can be cal- 
ibrated, then the results become much more inter- 
esting. Using contemporary priors, the limits range 
between tddm^ 33Hq^ and tddm > 43i7o"^ for DES, 



^ and tddm^ 66Hq ^ for our Deep survey, 



and tddm^ IOOHq^ and tddm^ 125Hq^ for our Wide 
survey. The variation between the lower values and 
higher values exhibits the range of possible constraints 
estimated using different nonlinear structure formation 



prescriptions. In all cases, the standard halo model gives 
the poorest constraint and the halo model modified to 
account for mass loss as the dark matter decays, as de- 
scribed in ? IIVI gives the most stringent constraint. The 
ability to exploit nonlinear power enables weak lensing 
to constrain unstable dark matter more stringently than 
contemporary methods using contemporary priors. 

Our most ambitious constraints come from assuming 
that Planck constraints on cosmological parameters other 
than r are available and that nonlinear structure forma- 
tion can be calibrated sufficiently to make full use of weak 
lensing data on nonlinear scales alongside Planck priors. 
These constraints are listed in the lower section of Ta- 
ble |T1 In this most ambitious scenario, the lensing con- 
straints on unstable dark matter are tddm^ 50Hq^ for 
DES, TDDM^ lOOi/o"^ for a Deep, JDEM/WFIRST-like 
survey, and tddm > ITOifp-^ for a Wide, LSST- or Euchd- 
like survey. Under these circumstances, weak lensing 
will provide the most stringent, model-independent con- 
straints on the decay lifetime of the dark matter particle. 

Numerous observational systematics need to be con- 
trolled for these instruments to achieve their statistical 
limitations and the theory on nonlinear structure growth 
in models with DDM must also be computed more rig- 
orously. However, our ambitious forecast limits exceed 
contemporary bounds on dark matter with invisible de- 
cay channels significantly, so that even if observational 
or theoretical systematics persist, weak lensing may yet 
provide the strongest limits on DDM. This opportunity 
(along with the related considerations in pro- 
vides a strong argument for a large-scale computational 
program to study the nonlinear evolution of density fluc- 
tuations in models with DDM. 

We conclude our results section, by addressing two re- 
maining outstanding issues. First, the most obvious stan- 
dard cosmological parameters that we may suspect to be 
degenerate with DDM are neutrino mass and the nor- 
malization of the matter power spectrum on large scales. 
Much of the constraining power of lensing comes from 
comparing contemporary potential fluctuations to those 
measured using the CMB at high redshift whereas signifi- 
cant neutrino mass also gives rise to scale-dependent sup- 
pression of the potential power spectra. Figure[2]provides 
anecdotal evidence that neutrino mass and DDM should 
not be so degenerate as to destroy constraining power be- 
cause they change the linearly-evolved potential spectra 
in different ways. Most importantly, the effect of DDM 
is strongly redshift dependent. Figure [5] shows projected 
confidence contours projected onto two-dimensional sub- 
spaces of our parameter space. Figure [5] shows that sig- 
nificant degeneracies do exist when only contemporary 
priors are used, and Planck priors suffice to break most 
of these degeneracies, leaving a slight degeneracy with 
neutrino mass as the most prominent. In the limit of a 
known neutrino mass (perhaps constrained by laboratory 
experiments), our constraints on DDM lifetimes improve 
by roughly ^ 25%, yielding a best-case constraint from 
a WIDE survey of tddm^ 210H^^. 
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FIG. 5; Marginalized, Icr constraint contours for dark mat- 
ter decay rate against other varying cosmological parameters 
in our model. In this illustrative example, we show contours 
from our linear power spectrum, ^max = 300, Wide survey 
calculations. The outer contours were computed using our 
contemporary priors while the inner contours were computed 
assuming Planck-level priors. This example is useful because 
degeneracies are most significant in the linear power spec- 
trum calculations. Notice that the degeneracy between dark 
matter lifetime and neutrino mass is relatively insignificant 
for contemporary priors, but is the most important remain- 
ing degeneracy with Planck prior constraints. The crosses 
designate fiducial parameter values. 



Finally, we have performed several direct searches 
through reduced cosmological and DDM parameter space 
to support our use of the Fisher matrix approach in the 
full parameter space. The need to reduce the param- 
eter space in direct searches is to limit computational 
effort. At present, with the nonlinear evolution in DDM 
models still uncertain, it does not seem fruitful to spend 
significant computational effort on forecasting. In a sim- 
ple, linear power model in which we scan the parameter 
space of r, cjm, Wf,, In(A'j^), and S^a, we find that our 
direct search decay rate bounds agree with the Fisher 
estimates to within ~ 20% for our Wide survey. More- 
over, the Fisher estimates are typically weaker than the 
projected direct search limits, largely because the Fisher 
matrix allows for degeneracies among parameters that 
would not be permitted on physical grounds. Our Fisher 
matrix forecasts are insensitive to the fiducial value of F, 
so long as the DDM lifetime is significantly larger than a 
Hubble time. 



VII. CONCLUSIONS 

We examined the utility of forthcoming, large-scale 
imaging surveys to constrain the lifetime of dark mat- 
ter decay into light daughter particles. Decaying dark 
matter can be disentangled from dark energy, because 
its primary observational signature is to reduce the po- 
tential fluctuation power spectrum, while dark energy 
is primarily constrained by geometric effects [s^, l87l |. 
DDM can be distinguished from massive neutrinos be- 
cause the suppression of power is a strong function of 
redshift. Assuming a null hypothesis of stable dark mat- 
ter, we found that utilizing only the information from 
linear scales 300) may suffice to place competitive 
limits on dark matter decay rates. Our linear-only, 
conservative forecast limits are F^ 13^0-^ for DES, 
F-i> 12H^^ for a Deep JDEM/WFIRST-like survey, 
and F-i> 18H^^ for a Wide LSST- or Euchd-like survey. 
The DES limit is slightly weaker than the best contem- 
porary constraints [l3, HO, H!! ; while the Wide limit 
is stronger. The constraint comes largely from a scale- 
and redshift-dependent reduction in convergence power. 
This indicates that without any further theoretical de- 
velopment, such instruments should provide competitive, 
complementary limits on unstable dark matter. 

In practice, lensing surveys will likely measure conver- 
gence spectra over a wide range of scales where nonlin- 
ear evolution is important to model {£ ~ 10^). The im- 
proved signal-to-noise from measurement on these scales 
should dramatically improve upon our conservative fore- 
casts. Assuming that we can use a halo model to map 
our linear spectra onto nonlinear spectra, we forecast lim- 
its that may be as stringent as F~^^ 50Hq^ for DES, 
F-i> lOOiJo-^ for a Deep JDEM/WFIRST-like survey, 
and F-i> 170iJo"^ for a Wide LSST- or Euchd-like sur- 
vey when combining with Planck satellite data. These 
forecasts are more restrictive than constraints available 
via other means. 

Of course, these more aggressive limit forecasts come 
with caveats. Developing surveys must control numer- 
ous nontrivial systematics to bring their lensing pro- 
grams to fruition and nonlinear evolution in models of 
DDM has not been explored well theoretically. Weak 
lensing systematics will be explored in support of the 
well-established goal of using lensing to constrain dark 
energy. A few groups have performed simulations with 
DDM models to understand its effect on dark matter ha- 
los and galaxy formation [l^ . However the number of 
samples and the simulated halo mass range are not yet 
sufficient to provide an adequate nonlinear fit to mat- 
ter power spectra. It is our hope that the opportunity 
to limit unstable dark matter using information on non- 
linear scales (see also will motivate researchers to 
explore nonlinear structure formation in models of un- 
stable dark matter. If this can be done, it may be pos- 
sible to limit the dark matter lifetime to be greater than 
hundreds of Hubble times. Furthermore, it would worth- 
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while to extend such calculations to models with small 
parent-daughter mass splittings. Such models introduce 
an additional parameter (the mass splitting) that deter- 
mines the recoil energies of the particles after decay. Such 
recoils may alter nonlinear structure at a level detectable 
with lensing data [3, ■ The future is promising for 
limiting the instability of dark matter using forthcom- 
ing, large-scale astronomical surveys. 
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